ValueFolder2='/sps/isc/rneveu/EtudeIRM/Nature/ResultsRetracted/';
VecResultsTemp=zeros(111,
for n3=0:11
    VecResults=zeros(2,3);
    for n1=1:161
        VecValueTemp=dir([ValueFolder2 'ResultsRetractedAuthors' num2str(n1) 'ArticleBefore.mat']);
        if size(VecValueTemp,1)==1
            load([ValueFolder2 VecValueTemp(1,1).name]);
            for n2=1:size(VecResultsRetracted,1)
                ValueTemp=find(VecResultsRetracted(n2,3:end-1)>0);
                ValueTemp=size(ValueTemp,2);
                if VecResultsRetracted(n2,2)>0
                    ValueTemp2=1;
                else
                    ValueTemp2=0;
                end
                if VecResultsRetracted(n2,end)>0
                    ValueTemp3=1;
                else
                    ValueTemp3=0;
                end
                if ValueTemp+ValueTemp2+ValueTemp3==0
                    VecResults(n2,1)=VecResults(n2,1)+1;
                else
                    if ValueTemp2+ValueTemp3==0
                        VecResults(n2,2)=VecResults(n2,2)+1;
                    else
                        VecResults(n2,3)=VecResults(n2,3)+1;
                    end
                end
            end
        end
    end
end

ValueFolderTemp='/sps/isc/rneveu/EtudeIRM/Nature/ResultsJournal/';
load([ValueFolder2 'AnalyseAuthorsRetractedArticle.mat']);
ValueFilter=0;
VecResults2=zeros(size(VecResults,1),3);
VecResults3=[];
VecResultsTemp=[1 2 3 4 5 10 0];
for n1=1:size(VecResultsTemp,2)
    ValueTime=VecResultsTemp(1,n1);
    for n2=1:size(VecAuthorsRetracted,1)
        ValueDate11=VecAuthorsRetracted{n2,3};
        if ValueDate11(1,1)>=1990
            ValueTemp=VecAuthorsRetracted{n2,2};
            load([ValueFolderTemp 'AnalyseMixedModelJournalRetractedTime' num2str(ValueTime) 'Filtered' num2str(ValueFilter) '.mat']);
            VecValueTemp=find(VecValueFigure(:,5)==ValueTemp);
            if size(VecValueTemp,1)>0
                VecValueTemp=VecValueFigure(VecValueTemp,:);
                ValueTemp2=find(VecValueTemp(:,4)==ValueDate11(1,1));
                if size(ValueTemp2,1)>0
                    VecResults2(n1,1)=VecResults2(n1,1)+VecValueTemp(ValueTemp2,2);
                    VecResults2(n1,3)=VecResults2(n1,3)+VecValueTemp(ValueTemp2,1);
                    VecResults2(n1,2)=VecResults2(n1,2)+VecValueTemp(ValueTemp2,3);
                else
                    if n1==1
                        VecResults3=[VecResults3;n2];
                    end
                end
            else
                if n1==1
                    VecResults3=[VecResults3;n2];
                end
            end
        else
            if n1==1
                VecResults3=[VecResults3;n2];
            end
        end
    end
end

VecResults4=zeros(7,3);
for n1=1:161
    ValueTemp=find(VecResults3==n1);
    if size(ValueTemp,1)==0
        VecValueTemp=dir([ValueFolder2 'ResultsRetractedAuthors' num2str(n1) 'ArticleBefore.mat']);
        if size(VecValueTemp,1)==1
            load([ValueFolder2 VecValueTemp(1,1).name]);
            for n2=1:size(VecResultsRetracted,1)
                ValueTemp=find(VecResultsRetracted(n2,3:end-1)>0);
                ValueTemp=size(ValueTemp,2);
                if VecResultsRetracted(n2,2)>0
                    ValueTemp2=1;
                else
                    ValueTemp2=0;
                end
                if VecResultsRetracted(n2,end)>0
                    ValueTemp3=1;
                else
                    ValueTemp3=0;
                end
                if ValueTemp+ValueTemp2+ValueTemp3==0
                    VecResults4(n2,1)=VecResults4(n2,1)+1;
                else
                    if ValueTemp2+ValueTemp3==0
                        VecResults4(n2,2)=VecResults4(n2,2)+1;
                    else
                        VecResults4(n2,3)=VecResults4(n2,3)+1;
                    end
                end
            end
        end
    end
end

VecResults5=[];
for n1=1:111
    VecValueTemp=dir([ValueFolder2 'ResultsRetractedAuthors' num2str(n1) 'ArticleBeforeAfter.mat']);
    if size(VecValueTemp,1)==1
        load([ValueFolder2 VecValueTemp(1,1).name]);
%         VecResults5=[VecResults5;VecResultsRetracted(:,2:end)'];
        VecResults5=[VecResults5;VecResultsRetracted(:,2)';VecResultsRetracted(:,end)'];
    end
end
VecValueTemp=find(VecResults5(:,1)<=15&VecResults5(:,2)<=15);
[ValueTemp,ValueTemp2]=ttest(VecResults5(VecValueTemp,1),VecResults5(VecValueTemp,2));
[ValueTemp3,ValueTemp4]=signrank(VecResults5(VecValueTemp,1),VecResults5(VecValueTemp,2));
VecResults6=[mean(VecResults5(VecValueTemp,2)-VecResults5(VecValueTemp,1)) sqrt(var(VecResults5(VecValueTemp,2)-VecResults5(VecValueTemp,1))) ValueTemp2 ValueTemp3]
                    
                
            
            
        